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We model coherent exciton transport in dendrimers by continuous-time quantum walks (CTQWs). For den- 
drimers up to the second generation the coherent transport shows perfect recurrences, when the initial excitation 
starts at the central node. For larger dendrimers, the recurrence ceases to be perfect, a fact which resembles 
results for discrete quantum carpets. Moreover, depending on the initial excitation site we find that the coherent 
transport to certain nodes of the dendrimer has a very low probability. When the initial excitation starts from 
the central node, the problem can be mapped onto a line which simplifies the computational effort. Further- 
more, the long time average of the quantum mechanical transition probabilities between pairs of nodes show 
characteristic patterns and allow to classify the nodes into clusters with identical limiting probabilities. For the 
(space) average of the quantum mechanical probability to be still or again at the initial site, we obtain, based on 
the Cauchy-Schwarz inequality, a simple lower bound which depends only on the eigenvalue spectrum of the 
Hamiltonian. 

PACS numbers: 71.35.-y, 36.20.-r, 36.20.Kd 



I. INTRODUCTION 

In recent years dendrimers have been an active field of 
research, both experimentally and theoretically, see, e.g. 
reference 1 or the Special Issue of Journal of Luminescence 
on the optical properties of dendrimers. 2 Dendrimers are hy- 
perbranched macromolecules with a very regular structure. 
Among a series of very interesting and crossdisciplinary ap- 
plications like drug delivery, dendrimers have been theoret- 
ically investigated as light harvesting antennae. 34,5,6 Apart 
from these theoretical works, there has been also a huge ex- 
perimental effort to probe transport processesii&2iI2iIi Here, 
dendrimers are synthesized in a self-similar fashion by hier- 
archically growing the dendrimer from a core^ Depending 
on the site of the light absorbing states, the transport of these 
excitations can be very efficient in some cases but also very in- 
efficient in others. At high temperatures the transport is inco- 
herent and can be described by a hopping process. 8 However, 
there is also experimental evidence for coherent interchro- 
mophore transport processes within the dendrimer. 1011 When 
modeling the transport, one usually assumes that the excita- 
tions are localized on the building blocks of the dendrimer, 
i.e., either on the chromophores or the segments connecting 
the chromophores. The latter leads to a different topology 
usually referred to as being a (finite) Husimi cactus^ 

There is a long-standing study of exciton transport in 
molecular aggregates, not only in polymer physics 14 but 
also in atomic^ and in solid state physics»A&Ai In solid 
state physics, a dendrimer of infinite generation is known as 
the Bethe lattice. The incoherent exciton transport in den- 
drimers can be efficiently modelled by random walks, see, for 
instancer&A2i2£. Here, the underlying topology of the den- 
drimer determines the dynamics of the exciton motion and the 
transport is described by a master (rate) equation. 

In this paper we will consider only the coherent transport, 
and we will model the dynamics by Schrodinger's equation. 
Interestingly, this is formally closely related to the master 
equation approach, where the transfer over the system is given 
by the connectivity matrix A of the dendrimerpiiS 2 * 2 ^! vide 



infra. This is also in close relation to Hiickel's (or LCAO, 
linear combination of atomic orbitals) theory, 2 - where the el- 
ements of the secular matrix are given by A. 



II. COHERENT EXCITON TRANSPORT ON GRAPHS 

We model the coherent transport of excitons on graphs by 
so-called continuous-time quantum walks (CTQWs) which 
are the quantum mechanical analog of continuous-time ran- 
dom walks (CTRWs). 

A graph is a collection of connected nodes and to ev- 
ery graph there exists a corresponding connectivity matrix 
A = (Aij), which is a discrete version of the Laplace op- 
erator. The non-diagonal elements equal —1 if nodes i 
and j are connected by a bond and otherwise. The diagonal 
elements An equal the number of bonds which exit from node 
i, i.e., An equals the functionality /j of the node i. 

Classically, a CTRW is governed by a master equation for 
the conditional probability, Pk,j{t), to find the CTRW at time 
t at node fc when starting at node j^&H The transfer matrix of 
the walk, T = (Tkj), is related to the connectivity matrix by 
T = — 7 A, where we assume for simplicity the transmission 
rate 7 of all bonds to be equal and we set 7 = 1. 

CTQWs are obtained by identifying the Hamiltonian of the 
system with the (classical) transfer operator (matrix), H = 
_iji ( 2i j 28 j 29 rpj ie w j 1G j e access ible Hilbert space is spanned 

by the basis vectors \j) associated with the nodes j of the 
graph. A state \j) starting at time to evolves in time as 
\j(t)) = V(t,t )\j), where U(t,t ) = exp(-iH(t - to)) 
is the quantum mechanical time evolution operator (we have 
set m = 1 and h = 1). 

The transition amplitude ak,j(t) from state \j) at time 
to state at time t reads then cekj(i) = (fc|U(t, 0)\j) and 
obeys Schrodinger's equation. Denoting the orthonormalized 
eigenstates of the Hamiltonian H = — T by \q n ), such that 
J2 n \Qn)(ln\ — 1> the quantum mechanical transition proba- 
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bility is 

ir kJ (t) = \a kij (t)\ 2 



-iA„i 



q n )(qn\j) 



(1) 



Note that classically ^ fc pk.j (t) = 1, whereas quantum me- 
chanically J^k \ a k.j(t)\ 2 — 1 holds. 



III. TOPOLOGY & CONNECTIVITY OF DENDRIMERS 

In the following we will consider dendrimers of generations 
G whose functionality is 3, i.e. all internal nodes of the den- 
drimer have 3 bonds, whereas the outermost nodes have only 
one bond. The stucture of such dendrimers is exemplified 
in Fig. [2 for dendrimers of generations G = 2 and G = 3. 
Note that the number of nodes belonging to the g-th genera- 
tion (g > 1) is 3 • 2 S_1 and that it grows exponentially with 
g. Moreover, the total number of nodes in the dendrimer of 
generation G is N = 3 • 2 G - 2. 





FIG. 1: Dendrimers of generation G = 2 (left) and G — 3 (right). 

The connectivity matrix of these dendrimers has a very sim- 
ple structure. One has An = 3 for all the nodes in the first 
G — 1 generations and An — 1 for the nodes in generation 
G. The bonds are represented by the off-diagonal matrix ele- 
ments Aij . Here, every node in generation 9 > 1 is connected 
to two consecutively numbered nodes in generation g + 1 and 
to one node in generation g — 1. 

The eigenmodes of such dendrimers were studied in^ 
There, for the dendrimers of generations G = 1 and G = 2, 
the eigenvalues and eigenvectors of A were explicitly calcu- 
lated. The eigenvectors determine the eigenmodes of the den- 
drimer, see e.g. Fig. 3 o&L. It was further shown that there 
are G + 1 nondegenerate eigenvalues, one of which is always 
Ao = 0. We used these results to check that our numerical 
diagonalization for the small dendrimers is correct. 



IV. TRANSPORT ON DENDRIMERS 

In the following section we present results for coherent 
exciton transport on dendrimers. The results were obtained 
by numerically determining the eigenvalues and eigenvectors 
of the corresponding connectivity matrix, using the standard 
software package MATLAB. 



A. Dynamics of an excitation starting at the center 

We start by focusing on an excitation which starts at the 
central node 1. For a G — 3 dendrimer, Fig.|2]shows snap- 
shots of the temporal development of transition probabilities. 
The nodes are numbered according to FigQ In this way the 
nodes 2 — 4 belong to generation g = 1, the nodes 5 — 10 to 
g = 2, and the nodes 11 — 22 to g = 3. Here, the excitation 
spreads over the whole dendrimer and already after a short 
period of time, t as 1.5, there is a considerable probability to 
be at the outermost nodes, which is larger than the probabili- 
ties to be at nodes belonging to other generations. The initial 
"speed" of the exciton is calculated based on the number of 
bonds (i.e. the chemical distance) between the initial node and 
one of the outermost nodes, divided by the time to reach that 
node. Here, we have 3/1.5 = 2. As we will show below, the 
dynamics of an excitation starting at the central node 1 over 
the dendrimer can be mapped onto a line. For an infinite line, 
the transition probabilities can be expressed by Bessel func- 
tions of the first kind. 22 From this the speed can be calculated; 
it is at all times equal to 2. For finite lines, the initial speed is 
still equal to 2. Thus our results here are in accordance with 
what can be expected for waves propagating through a regular 
graph. 
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FIG. 2: Snapshots of the transition probability for a G = 3 den- 
drimer at different times t. The initial excitation is at the central 
node. 
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Remarkably, for the G — 2 dendrimer the transition prob- 
abilities are fully periodic when the coherent excitation starts 
from the central node 1 (the same holds for the G — 1 den- 
drimer, too). In this case we obtain, based on the analytically 
determined eigenvalues and eigenvectors that all iTk.i{t) have 
the following, periodic form: 



Hgig and are given by 



7Tfc,l(t) 



E« 



(k\q n ){q n \±) 



= a { k) + a[ k) cos(2t) + a { 2 k) cos(3i) + a {k) cos(5t).(2) 

Note that due to rotational symmetry, the transition probabili- 
ties from the central node to nodes belonging to the same gen- 
eration are equal. Because of this we only have to list three 
different transition probabilities. We hence choose, exemplar- 
ily, the nodes 1, 2, and 10, for which we find the coefficients 
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If follows that there is a perfect revival of the initial state 
for each t = 2mr where n € N. At t = (2n — l)ir the main 
part of the excitation is equally distributed among the nodes of 
the outermost generation. This revival of the initial probabil- 
ity distribution resembles results obained for continuous 31 ' 32 
and discrete quantum carpets. 22 - 33 For discrete quantum car- 
pets, the revival is only perfect for small cycles of length 
N = 1,2, 3, 4, and 6. When the cycles become larger, there 
are only partial revivals of the initial state. The same can be 
observed for dendrimers of generation G > 3, as will become 
clear in the following. 



B. Mapping onto a line 

Due to the rotational symmetry, the dynamical problem of a 
dendrimer initially excited at its center 1 can be mapped onto 
a linear segment. In the classical, CTRW case, this result has 
been widely used, see Ref. 5 . For a CTQW on a graph consist- 
ing of two Cayley trees this has been done in, e.g.^i^. For 
each generation with g > 1 we introduce linear combinations 
of states 



Is) 



=y>>- 

oa-i ^ 1 ' 



(3) 



keg 



The state of generation 0, \go), is identical to that of the central 
node 1 1) . A CTQW through the generations is now governed 
by a new Hamiltonian H which is defined by 



H gg , = (g\U\g') 



(4) 



Given H and the construction of the generation states \g), H 
is a real and tridiagonal matrix whose elements obey H gg i — 



H, 



g,g±i 



(s|H| ff ) = f g , 
(3|H| 5 ±1> = 



(5) 



3 • 29~ 

kegjeg±l 

V3 for H i = H w 
a/2 else. 



(6) 



Here, f g is the functionality of the nodes in generation g and 
all other elements of H are zero. In this way the original 
eigenvalue problem which grows exponentially with G has 
been reduced to an eigenvalue problem which grows only lin- 
early with G. The G + 1 eigenvalues of the new matrix H 
are the nondegenerate eigenvalues of the original matrix H. 
In terms of eigenmodes, these are the modes of the full den- 
drimer, where, in a classical picture, 30 nodes belonging to the 
same generation move in the same direction. 

Now, the transition amplitude between different genera- 
tions can be written as (g\ exp(— iHt)\go) . It is easy to show 
that the transition probabilities for a line of 3 nodes obtained 
from the G = 2 dendrimer are, see Eq.0, 



jcx P HHi)| 5o )| 2 =^7T M (t). 

keg 



(7) 



The cumulative transition probabilities for the G = 2 and 
the G — 3 dendrimer to go from one generation of the den- 
drimer to another one when starting at g — are shown in 
Fig- E] For generations G > 3, we do not find perfect revivals 
of the initial probability anymore. 




FIG. 3: (Color online). Cumulative transition probabilities that a 
coherent excitation starting at the central node is at time t at a node of 
generation g. We present the situation for dendrimers of generations 
(a) G = 2 and (b) G = 3. Note the complete revival at t = 2n for 
the G = 2 dendrimer. 
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C. Dynamics of excitations starting off-center 



branches of the dendrimer stays very low. 



If the initial excitation is placed at one of the nodes of the 
outermost generation g = G of the dendrimer, the picture 
changes. Snapshots of the transition probability for a G — 
3 graph where the excitation starts at the outermost node 22 
are shown in Fig.|4] Even classically the propagation through 
the dendrimer gets to be much slower than in the previous 
case, see e.gpi&^i^. Nevertheless, eventually the excitation 
will classically propagate through the whole graph and in the 
long time limit the probability will be equipartitioned among 
all nodes. 
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FIG. 4: Snapshots of the transition probability for a G = 3 den- 
drimer at different times t. The excitation starts at the outermost 
node 22. 

Quantum mechanically this effect is even more dramatic. 
For the times shown in Fig. |4] the main fraction of itk,j(t) 
stays in a small region closely connected by bonds to the ini- 
tial node j, and the transfer to other sites is highly unlikely. 
For the G — 3 dendrimer shown in Fig. [2 the excitation is 
located with very high probability on the initial node 22 and 
on the nodes 21 and 10. On the time scales shown in Fig.|4] 
the probability for nodes outside the branch (consisting of the 
nodes 1, 4, 9, 10, 19, 20, 21, and 22) to be excited is very low . 

As we show in the following section, also at long times 
the limiting probability for the excitation to reach the other 



D. Long time limit 

Quantum mechanically the temporal development is sym- 
metric to inversion; this prevents Ttk,j{t) from having a defi- 
nite limit for t — > oo. In order to compare the classical long 
time probability with the quantum mechanical one, one usu- 
ally uses the limiting probability (LP) distribution 36 



Xk,j = T li] ^^ / lU '^/('i 



(8) 



which can be rewritten by using the orthonormalized eigen- 
states of the Hamiltonian, \q n ), 



Xk,j 



= lim 



dt 



£<*!« 



-Oft 



Qn)(q n \j) 



A„.A„ 



\qn)(qn\j)(j\q m )(q m \k). 



(9) 



Some eigenvalues of H might be degenerate, so that the sum 
in Eq. (|9jl can contain terms belonging to different eigenstates 
\q n ) and \q m ). 
We can use the Cauchy-Schwarz inequality 37 



dt f(t)g(t) 




dt \f(t)\ 




dt\g(t)f 



(10) 

to obtain a lower bound for the LP. With f(t) = ak,j(t) and 
g(t) = 1, the time integral in Eq. l|8} fulfills the inequality 



dt\a k ,j(t)\ 2 > 



T 



dt a k ,j(t) 



(11) 



This results in 



Xfe, 



> lim 



-iXnti 



q n ){q n \i) 



(12) 



The only term in the sum over n in Eq. (1121 which survives 
after integration and taking the limit T — > oo is the one with 
Ao = 0. In terms of eigenmodes, the eigenstate correspond- 
ing to Ao is the one for which in a classical picture the den- 
drimer moves as a whole in one direction. The corresponding 
eigenvector can be written as |tjo) = l/v^"X]^Li Since 
the states \j) form a complete, ortho-normalized basis set, i.e. 

(k\j) — Skj, we get (fc|go) = (qo\j) = l/VN and therefore, 
with Eq. O, 



2 1 

Xk,j > \ (k\q )(qo\j)\ = j^- 



(13) 



Classically, the LP is equipartitioned among all the nodes, 
i.e. lirrit-foo Pk,j (t) — 1/iV for all nodes. Quantum me- 
chanically this is not the case. However, for regular graphs 
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like rings and for initial conditions localized on one node, 
the quantum mechanical LPs are almost equipartitioned, i.e., 
there are only a few nodes whose LPs differ from the rest. 24 
For more complex structures, the LPs display large variations, 
depending on the nodes, as we proceed to demonstrate. 
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FIG. 5: Limiting probabilities for dendrimers with (a) G = 3 and 
(b) G = 4. The white lines indicate the limiting distributions for the 
dendrimers of generations smaller than G. 

Figure |5] shows the LPs Xk,j as a contour plot where the 
axes are labeled by the nodes k = 1, • ■ • ,N and j — 
1, • • • , N, Bright shadings correspond to high values of the 
LP, whereas dark shadings correspond to low LPs. The diag- 
onal has high values, meaning that excitations starting at any 
node j have a high LP to be at node j again. Note that all 
values of Xk,j are clearly larger than the lower bound given 
in Eq. Jl 31 . which is (1/22) 2 for the G = 3 dendrimer and 
(1/46) 2 for the G = 4 dendrimer. 

The structures of the LP distributions of dendrimers are 
self-similar generation after generation. The white lines in 
Fig.[5]are a guide to the eye, indicating this fact. 

Furthermore, different nodes k and I may have the same LP, 
Xk.j = Xi.j- We hence combine all the nodes having (up to 
our numerical precision, 10~ 10 ) the same LP into a cluster. 
Note, however, that the separation of the nodes into clusters 



depends on the initially excited node, namely on j. For an ex- 
citation starting at the center (node 1), the clusters correspond 
exactly to the different generations of the dendrimer. In the 
general case, when starting from a non-central node, we still 
find from Fig.[3]that nodes belonging to the same cluster also 
belong to the same generation (the converse is not necessarily 
true). 




FIG. 6: (Color online). Upper two panels: Limiting probabilities for 
the G — 2 (top row) and the G = 3 (bottom row) dendrimer when 
the excitation starts from the central node 1 [(a) and (b)] and from 
one of the outermost nodes [(c) from node 10 and (d) from node 22]. 
Lower panel: Groups of nodes (clusters) having the same limiting 
probability \k,j for a dendrimer of generation G — 3. In (e) the ex- 
citation starts from node 1 and in (f) from node 22. Nodes connected 
by thick (red) lines belong to the same cluster. 

We now focus on the two special cases of initial excitations 
given in the previous sections. Figure [6]presents the LP dis- 
tributions for the G = 2 and G = 3 dendrimers; in both cases 
the excitation starts either at the central node 1 or at one of the 
outermost nodes (node 10 for G = 2 and node 22 for G = 3). 

Having the initial excitation at the central node results in a 
totally symmetric LP distribution, i.e. in each generation all 
nodes have the same LP. In this case we can again simplify 
the problem by considering only the cumulative probabilities. 
The LP given in Eq. simplifies, because the cumulative 
probabilities contain only nondegenerate eigenvalues, i.e. 



Xg,g' 



(14) 



Here, the \q n ) are the eigenstates of the reduced Hamiltonian 
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H, which does not have degenerate eigenvalues. For the to- 
tally recurrent case of an excitation starting at the center of the 
G = 2 dendrimer, the LP distribution follows directly from 
Eq. (0 as 



Xk,i 



,0) 



(15) 



Note that for an excitation starting at the center the number of 
different clusters equals G + 1. 

For an excitation starting from one of the outermost nodes, 
the LPs are less regular. However, also in this case the LPs 
cluster. Figure |6jf) shows the situation for the G = 3 den- 
drimer when starting from node j = 22. Here, nodes belong- 
ing to the same cluster are connected by thick (red) lines. The 
cluster structure beyond the central node 1 can be understood 
on simple symmetry grounds. The structure on the side of the 
initial node, however, is more complex. Here, nodes 22 and 21 
have the same LP and form one cluster. The same holds (re- 
markably) for nodes 9 and 10 (another cluster) and for nodes 
19 and 20 (yet another cluster). 
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FIG. 7: (Color online). Clusters of the same limiting probability 
Xk,j in the branch with the initial excitation for a dendrimer of gen- 
eration G. Nodes connected by thick (red) lines belong to the same 
cluster. 

For larger dendrimers, while the general cluster pattern is 
preserved, some details change. Figure shows the situation 
for a dendrimer of dimension G = 5, for an excitation starting 
at a peripheral node. Again we indicate clusters by connecting 
nodes by thick (red) lines. A change to be noticed is that for 
G > 5 the initially excited node does not form anymore a 
cluster with its next-nearest node of the same generation. It 
appears as if such two nodes belong to the same cluster only 
when the dendrimer has G < 4. Thus the total number of 
clusters is N c = (G 2 + G + 6) /2 for G > 5 and N c - 1 
for G < 4. Moreover, the study of the dendrimer with G = 6 
shows an analogous situation to our findings for G = 5. 



E. Averaged probabilities 



are periodic. However, for larger dendrimers and/or different 
initial conditions this does not have to be the case. 

Classically, there exists a simple expression for the average 
probability to be still or again at the initially excited node. 
One has (see e.g. Refi^Sl): 

N 1 n 

= n X>wM = N ^2 CXP ( ~ A "*)' (16) 
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n=l 



This result is quite remarkable, since it depends only on the 
eigenvalue spectrum of the connectivity matrix A but not on 
the eigenvectors. 
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FIG. 8: (Color online). Comparison of the averaged probability to 
be still or again at the initial node with the lower bound given in 
Eq. j2Ql for generations (a) G = 2, (b) G = 3, and (c) G = 4. 

Quantum mechanically, the average is given by 



1 N 



(17) 



By using the Cauchy-Schwarz inequality we obtain a lower 
bound for 7f(i), 



One interesting question to ask is, what is the probability, 
TTj,j(t), to be at the initial site after some time tl As shown 
above, for the G = 2 dendrimer 7Tfc i(i) and hence 7ri,i(t) 



N N 



i=i 



(18) 
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With Eq. Q we have 

E a ^ = EE e " itAB <«"iJ>oifc> = E e_itAB - (19) 

i i n n 

and therefore 

^ 772 E ex p [ - *( A « - A ™)^ ■ < 2 °) 

In analogy to the classical case, the lower bound in Eq. i20\ 
depends only on the eigenvalues and not on the eigenvectors 
of A. Note that for a CTQW on a simple regular network 
with periodic boundary conditions the lower bound in Eq. i20i 
becomes exact. 

Figure [8] shows the exact value of the average probability 
W(t) as well as the lower bound given in Eq. J20i for den- 
drimers of three different generations. One notes immediately 
that the fluctuations of the lower bound are larger that the ones 
of the exact value. However, the qualitative agreement be- 
tween the two curves is remarkable: The positions in time of 
the extrema of the lower bound and of the exact curve almost 
coincide. Especially the strong maxima of W(t) are well re- 
produced by the lower bound. 

In fact, fluorescence spectroscopy experiments like, e.g., 
the ones of Varnavski et al*&lL allow, in principle, to distin- 
guish whether the transport is rather classical or rather quan- 
tum mechanical. In order to identify excitonic coherence, 
Lupton et al. compared a coupled harmonic oscillator model 
to photoluminescence spectra. 40 Another suggestion is the one 
by Heijs et al, in which pump-probe and fluorescence exper- 
iments determine the efficiency of the excitation transfer, 18 
from which one can infer the underlying transport mechanism. 
We close by noting that such models are, as a rule, quite ideal- 



ized, so that one must be careful in relating the models to the 
experimental findings. 



V. CONCLUSIONS 

In conclusion, we have modelled the coherent exciton trans- 
port by continuous-time quantum walks. The transport is ex- 
clusively determined by the topology of the dendrimer, i.e. by 
its connectivity matrix. For the G = 2 dendrimer the transport 
is completely recurrent when the central node is initially ex- 
cited. In this case the quantum mechanical transition probabil- 
ities are fully periodic. For larger dendrimers and/or different 
initial conditions we observe only partial recurrences. 

To compare these results to those of the classical (incoher- 
ent) case, we calculated the long time average of the transition 
probabilities. Depending on the initial conditions, these show 
characteristic patterns. Furthermore, the limiting probability 
distributions can be characterized by clusters of nodes having 
the same limiting probabilities. 

Furthermore, we calculated a lower bound for the average 
probability to be still or again at the initial node after some 
time t. This lower bound depends only on the eigenvalues of 
A and agrees qualitatively well with the exact value; espe- 
cially the maxima of lf(t) are well reproduced by the lower 
bound. 
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